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Abstract 

Zipoy used spheroidal coordinates to construct a family of static 
axisymmetric gravitational fields included in Weyl's class. We calcu- 
late the mass moment source for the dipole solution using techniques 
previously developed for magnetic and angular momentum densities 
and compare it with the Newtonian analog studied by Bonnor and 
Sackfield. 
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1 Introduction 



In [1] Bonnor and Sackfield applied potential theory to analyse some 
metrics of the family obtained by Zipoy [2]. They calculated the 
mass density for the monopole solution from the discontinuities of the 
derivatives of the Weyl potential [3] , taken as a Newtonian potential, 
and also they calculated the moment density for the dipole solution 
from the discontinuity of the potential itself, considering these jumps 
as surface sources for the gravitational fields. While their result for 
the monopole solution was also grounded on the thin surface layer for- 
malism presented by Israel in [4] and therefore it had an interpretation 
within the framework of General Relativity, their analysis of the dipole 
solution was carried out using only a Newtonian flat spacetime analog, 
concluding that a sheet of dipoles was the only source present. In this 
paper a justification for that result will be given using a modification 
suited for mass dipoles of the formalism derived previously in [5] and 
[6] for rotation and electromagnetic moments in General Relativity. 

Section 2 is devoted to the Weyl static axisymmetric family of vac- 
uum solutions of the Einstein equations using in its derivation the dif- 
ferential form approach applied in [7] and [8] to vacuum fields. Zipoy's 
spheroidal metrics are dealt with in Section 3. In Section 4 use is made 
of the formalism introduced previously to provide a way of calculating 
dipole densities out of the discontinuities of the metric functions. This 
method is applied to the Zipoy dipole solution in Section 5. A brief 
discussion of the results is provided in Section 6. 

2 The Weyl Static Family 

The Weyl static axisymmetric vacuum solutions were discovered soon 
after the publication of the General Theory of Relativity [3] . It is the 
general solution for metrics with this symmetry and allows the con- 
struction of metrics from solutions of a flat spacetime Laplace equa- 
tion. 

We shall introduce the relevant equations for the problem in the 
form of a exterior system. Our start point is an orthonormal vierbein 
{6'^ , 9^ , 6'^ , 0^} such that the metric takes this form: 

4^ = -e^ + 9^ + 0"^ (^e"^ + (i) 
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The forms 6^ and 6^ lie on the space spanned by the orbits of the 
two Kihing fields, dt and d^, while the other two forms, 6^ and 9^, lie 
on the orthogonal space. 

The equations stating that the frame is torsion- free (2-5), together 
with its integrability conditions (6-7) and Einstein's equations (8-9) 
can be shown to be [7] , [8] : 

de^ = aA (2) 
d0^ = {b-a)Ae^ (3) 

de^ = -uAe^ (4) 
de^ = uA (5) 

db = (6) 
da = (7) 
d*b + bA*b = (8) 

d * a + b A *a = (9) 

The symbol * stands for the Hodge duality operator in the 6'^ — 6"^ 
space, that is = 6^ and = —9^. The equations governing the 
connection u in the 0^ — 9^ space have not been written since they can 
be integrated by quadratures after the previous exterior system has 
been solved and are not of interest for our purposes here. 

The equations 2-8 can be formally integrated to yield the metric 
in canonical coordinates {t,(p, z, p} and therefore the meaning of the 
one-form b as the differential of the logarithm of the volume of the 
9^ — 9^ space is clear: 

a = dU b = d\np *b = -p''^dz (10) 

S = -e^^dt^ + e-2^{p2#2 + e'^'^idp'' + dz'')} (11) 

The remaining equation, 9, is analagous to the Raychaudhuri equa- 
tion for the accelerarion of a fluid. Therefore a can be viewed as an 
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acceleration one-form. After substitution of 10 it takes the form of a 
reduced three dimensional Laplace equation: 

d*dU + p-^dp A*dU = (12) 
Or in explicit Weyl canonical coordinates: 

Upp + U,, + ^Up = (13) 

As we had already stated, Newtonian potential functions, U, yield 
after substitution in the metric 11 exact solutions of Einstein's equa- 
tions. But the analogy finishes here and no further interpretation can 
be brought from this fact. As an example, the classical potential for 
a finite rod with constant density produces the Schwarzschild metric, 
the relativistic monopole solution. 

If we want the metric to be asymptotically flat, then U should be 
at most monopolar at infinity, that is: 

We shall restrict to these metrics from now on. 



3 Zipoy's Metrics 

Zipoy's metrics belong to Weyl class [2] and are constructed from so- 
lutions of the Laplace equation in oblate spheroidal coordinates {r, 9} 
instead of using pseudocilindrical coordinates {p,z}: 

z = rcos9 p = Vr'^ + a'^sinO (15) 
where a is a constant and the new coordinates range as usual: 

0<r<oo 0<6l<7r (16) 

Note that events with coordinates (t, (p,r = 0, 9) and (t, (p,r = 
0,-jT — 9) have the same canonical coordinates. Therefore we shall 
identify them instead of attempting more complicated interpretations. 

In Bonnor and Sackfield's paper, is a latitude angle instead of a 
colatitude, as it is taken in this paper, and also we make use of the 
radius r instead of the coordinate u related to the previous ones by r = 
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asinhu. We have done this since their interpretation as asymptotic 
spherical coordinates is simpler. 

The reduced Laplace equation in these coordinates takes the fol- 
lowing form: 

(r^ + a'^)Urr + ^rUr + Uee + cot OUe = (17) 

We can construct solutions for this equation of the form Un = 
Rn{r)Pn{cos6) where P„ stands for the Legendre polynomials. The 
index n runs from zero to infinity allowing us to classify the Newtonian 
potentials U and the metrics derived from them using a multipole 
expansion. We shall be concerned with the solution for n = 1, Zipoy's 
dipole solution. 



4 Dipole Sources 



Since the aim of this paper is just the interpretation of Zipoy's dipole 
solution, we shall restrict ourselves to asymptotically dipolar static 
axisymmetric vacuum solutions of the Einstein equations. For this we 
mean that the metric function U must be: 

Mcos6i ^ / 1 \ 
u = —:^ + o (18) 



^2 V 

where M is the total mass dipole. 

Now, we shall proceed as in [5] and [6] and consider the two formal 
expressions for the acceleration form arising from the integration of 
equations 7 and 9: 

a = dU = -p-^ * dA (19) 

The asymptotic behaviour of U imposes a condition at infinity for 
the newly introduced function A: 

^ = ^ + (20) 

In order to have a dipole source for the gravitational field, it is 
assumed that the function U is discontinuous across a closed surface 
S. This condition on S means no restriction, since the difference 
between the values of U on either side of the surface can be null on 
some parts of it. 
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Taking into account both expressions for the acceleration, we inte- 
grate the scalar product of a and the differential of the Weyl coordinate 
z with a factor e^. The domain of integration will be the space V3 

orthogonal to 9^, whose metric is g g + 9^ ^ 9^. The square root 
of the determinant of this metric is ^yg = pe^^~^^ . 

= / Vge^ <[a + *{*a)],dz > dx^dx^dx^ = 
JV3 

= I Vae^ < [dU + p'^ * dA],dz > dx^dx^dx^ (21) 

The Wcyl coordinate z satisfies the following differential equation 
arising from 8, since ~ ^ \f9'- 

= (22) 

Therefore the integrand in 21 can be shown to be a total derivative 
and the expression can be reduced to a surface integral. 



= / ^^,{^U^ g^''d^z^\nv\K.d^z)dx^dx^dx^ (23) 

Actually, since U is discontinuous, the integral has to be split into 
two pieces, and ^ respectively the outer and inner sections of 
^3 referred to the surface S. The boundary of is formed by the 
surface S and the boundary of consists of the same surface and the 
2-sphcrc at infinity S'^(oo). The integral over 5*^(00) can be straight- 
forwardly carried out, taking into account the asymptotic conditions 
18 and 20 and that the metric is assumed to be asymptotically flat: 

Hence the integral defining the total mass dipole is just: 

47rM = / dS^We^ g^^'riad^zY - [ dS' {Ue^ g'"'n.d^z}~ (24) 
Js+ Js- 

Due to the fact that the metric is discontinuous -the potential U 
is also part of the metric and of the unitary normal, n, to the surface 
S-, it seems in principle not possible to reduce the previous equation 
to one integral of the jump of U. But it can be shown, for instance, 
that, if the discontinuity takes places on a disk on the plane z = 0, 
metric factors like cancel each other and therefore the differential 
element of mass dipole is perfectly defined. 
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5 The Dipole Solution's Source 



The aim of this paper was the interpretation of Zipoy's dipole metric 
with the formalism developed in the previous section. As it has al- 
ready been said, this metric arises from the substitution of the dipolar 
solution of the Laplace equation in oblate spheroidal coordinates in 
the Weyl metric. This solution can be shown to be: 

U = -f{l-- arctan(-)} cos 9 (25) 
a r 

where 7 is a constant. 

It is easily seen that this function behaves for large values of r like 
a dipolar potential whose moment is M = 570^. 



1 2 cos 9 
-7a , - , 

The remaining metric function, although it will not be needed is: 



t^ = >^^ + 4) (26) 



^ j ■ — ' - — - 1 exp{-i7^ tan^ 9 - j'^ sin^ e[arctan(-)]^} 



+ cos^ 9 ' 

r"^ + ] ^ ' ^ V 



(27) 

As we have identified on the surface of null r coordinate points 
whose colatitude is 9 with those that have ir — 9, it happens that, 
when we reach that surface from the subspace z > 0, the potential 
takes the value U = ^cos9, but if we arrive with 2; < 0, then it is 
equal to U = —jcosd, 9 ranging from to | after the identification 
of events [1]. 

Therefore the potential meets a jump at the disk z = 0, p < a 
(r = O,0<f): 



[U] = 27 cos 9 = 2j^l-^ (28) 
The disk z = has the following surface element and normal: 

dS = pe^-'^^ dpdcj) n = e^'^'d^ (29) 

Therefore the product eP n^d^zdS does not depend on U and takes 
the same value on both sides of the disk S, as we anticipated in the 
previous section. 
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Hence the integral defining the total mass dipole can be written 

as: 




dS \U] e' 



47r Jo 




(30) 



The mass dipole takes the value expected from the asymptotic 
form of the metric and the differential element of the mass moment is 



in full accordance with Bonnor and Sackfield's Newtonian result. 



In the previous section the differential surface element for the mass 
moment source of Zipoy's dipole metric is constructed from the dis- 
continuity of the metric function U. This result coincides with the one 
obtained by Bonnor and Sackfield [1] using classical potential theory, 
since the curved spacetime elements in the integral 30 arising from the 
metric cancel each other. Therefore the full rclativistic theory and its 
classical Newtonian aproximation yield the same output and Bonnor 
and Sackfield's result is confirmed within the framework of General 
Relativity. 

Something similar happens with the angular momentum density 
for Kerr's metric [5], where oblate spheroidal coordinates are also em- 
ployed. As it is stated in [1], the potential U is due only to the sheet of 
dipoles described in 31, that is, there is no contribution of higher order 
multipole sheets to the gravitational field and therefore the source is 
fully described, at least in the Newtonian limit. Again this is the case 
for Kerr's metric. 

It remains to be seen if this also happens when the full relativistic 
theory is taken into account. 
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6 Discussion 
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